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ABSTRACT 

The dominant source of electromagnetic energy in the Universe today (over ultraviolet, optical and 
near-infrared wavelengths) is starlight. However, quantifying the amount of starlight produced has 
proven difficult due to interstellar dust grains which attenuate some unknown fraction of the light. 
Combining a recently calibrated galactic dust model with observations of 10,000 nearby galaxies we 
find that (integrated over all galaxy types and orientations) only (11 ± 2)% of the 0.1 fim photons 
escape their host galaxies; this value rises linearly (with log A) to (87 ± 3)% at 2.1 fim. We deduce 
that the energy output from stars in the nearby Universe is (1.6 ± 0.2) x 10 35 W Mpc~ 3 of which 
(0.9±0.1) x 10 35 W Mpc -3 escapes directly into the inter-galactic medium. Some further ramifications 
of dust attenuation are discussed, and equations that correct individual galaxy flux measurements for 
its effect are provided. 

Subject headings: galaxies: spiral - galaxies: structure - galaxies: photometry - galaxies: fundamental 
parameters - ISM: dust, extinction 
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1. INTRODUCTION 

The cosmic spectral energy distribution (CSED; e.g., 
Primack, Bullock & Sommerville 2005) provides a de- 
scription of the current total (electromagnetic) energy 
output of the Universe over all wavelengths. In the 
ultraviolet, optical and near-infrared wavebands the 
CSED is dominated by starlight and its measurement 
can be used to constrain the current stellar mass den- 
sity and cosmic star-forma tion rate as well as models 
of galaxy formation (e.g ., iBaldry fc GlazebrooU 120031 : 
iHopkins fc Beacoml 120061 ). The CSED is measured 
by constructing the galaxy luminosity function (GLF) 
dSchechterl Il976l : iFeltenl [19771 IDriver et alJ l2005h at a 
specified wavelength (or bandpass). The first moment of 
the GLF (extrapolated to bright and faint magnitudes), 
gives the total luminosity density at this wavelength and 
hence provides a single datum on the CSED. 

Constructing the full CSED therefore requires accu- 
rate measurements of the GLF in a variety of band- 
passes. However, galaxies contain dust, whic h, while neg- 
ligible in terms of mass (jDriver et al.l 120071 ) . attenuates 
some unknown fraction of the starlight before i t exits 
a galaxy into the inter-gala ctic medium (IGM) (|Seared 
119311 iGiovanelli et al.lll995| ). The severity of this ef- 
fect has proven difficult to quantify (jDisnev et al.l [19891 
IValentiin| [l990: Bu rstein et al.|[l991l ). leading to large un- 
certainties in individual galaxy flux measurements and 
consequently a systematic underestimation of the lumi- 
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nosity density (or individual CSED measurements). The 
degree to which the CSED is un derestimated will, of 
course, be w avelength dependent (jCardelli et al.l 119891 : 
Calzctti 2001) and it will depend critically on the amount 
and distribution of the interstellar dust grains within the 
host galaxy (and an individual galaxy's orientation to our 
line-of-sight). 

Since the heated debate of the 90s, direct evidence 
from a number of methods have led to the perspective 
that galaxies are predominantly optically thin (at least 
in their outer region s). In particular, evidence from 
overlapping galaxies (j White et al.l 12000 ). have strongly 
supported the stance that gala xies are optically thin , 
at least in the inter-arm regions (|Holwerda et al J 120071) . 
Further detailed modeling of extensive optical data on 
edge-on galaxies also appeared to support the view that 
galax ies were optically thin throughout (jXilouris et al.l 
1999). However, the optically thin case has proven dif- 
ficult to reconcile with the observed high level of far- 
infrared emission which is presume d to arise from dust 
reradiating the attenuated starlight jlBianchi et al.l 120001 : 
iPopescu et all 120001 : iMisiriotis et alJl200lT >. The resolu- 
tion to this conflict may lie in a more complex dust dis- 
tribution whereby galaxies may contain both optically 
thick (core and arm ) regions an d optic ally thin (inter-arm 
and ou ter) regions. I Tuffs et ail ()2004l ) and lPopescu et all 
(2000) (hereafter TP) have advocated a three component 
dust model (optically thick inner disc, thin outer disc and 
clumpy components), which is capable of reproducing the 
detailed multiwavelength surface photometry from UV to 
far-IR of edge-on gal axies such as NGC891 and the other 
galaxies mo delled bvlXilouris et al.l (|1999f h 

Recently, IDriver et all (|2007l ) identified evidence for 
strong and inclination-dependent attenuation in a large 
sample of discs and bulges which is not anticipated 
(or reproduceable) in models with a purely optically 
thin dust distribution (e.g., see fig 5 in Popescu & 
Tuffs 2007). Similar and related results, albeit lack- 
ing bulge-disc decompositions, have also now been 
reported for SPSS data ( e.g., Choi, Park fc Voge- 
ley 2007; iShao et all 120071 : lUnterborm fc Rvdenl 120081: 
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iMaller et all l2008t iPadilla fc Straussl I20M) . Using the 
TP-model we were able to reproduce the attenuation- 
inclination relation for both discs and bulges and to con- 
strain the mean central _B-band face-on opacity (the only 
free parameter) for the galaxy population at large: Tg = 
3.8 ± 0.7. In this Letter we explore the implications of 
this result on estimates of the CSED and ask whether the 
high value for the central opacity can be reconciled with 
the far-IR output from the nearby galaxy population. 
Throughout we adopt a standard flat cosmology with 
n M = 0.3, n A = 0.7 and H = 70 h 70 km s" 1 Mpc" 1 . 

2. DUST AND ITS IMPACT ON GALAXY PHOTOMETRY 

The dust model we adopt here is described fully in 
a sequence of papers (IPopescu et al.l 120001 : iTuffs et"aL1 
20041 iMollenhoff et all 120061 ) andis summarised by 
Popescu & Tuffd (|2007t l. In brief the model incorpo- 
rates three distinct components: an extended optically 
thin dust disc associated with the neutral hydrogen and 
older stellar population (i.e., the outer disc and inter- 
arm regions); a less extended optically thick dust layer 
in the spiral arms associated with the molecular hydro- 
gen and younger stellar population (i.e., primarily the 
inner disc); and a clumpy component (representing star- 
forming molecular clouds). Ideally the optically thick 
component should be distributed according to a spiral 
pattern, to better mimick the observed variation in opac- 
ity between arm and inter-arm regions. For the purposes 
of this work this distinction, between an additional disc 
of uniform opacity or one with a built in spiral density 
pattern is not particularly relevant as both imply a uni- 
form high opacity disc in the central bulge-dominated 
regions. 

In iDriver et a l. (2007) we constrained the TP-model's 
single free parameter, the central face on opacity in £>, 
which was found to be = 3.8 ± 0.7. Briefly, this con- 
straint was obtained by measuring the luminosity func- 
tions of galaxy bulges and discs at various inclinations 
and comparing the dependence of the turn-over point on 
inclination (i.e., th e M*-cos(t) relatio n) to predictions of 
the TP-model (see IDriver et al.ll2007l for full details). 

Having constrained the TP-model we can now deter- 
mine the inclination-dependent attenuation correction at 
any wavelength. From Fig. 1 we see that the attenuation 
for bulge starlight can be as high as 2 B mag (i.e., only 
6% of the photons escape), and as high as 1.2 B mag 
for discs (i.e., only 33% of the photons escape), depend- 
ing on wavelength and inclination. While the level of 
attenuation of the disc is relatively consistent with pre- 
vious estimates, the bulge attenuation is a surprise and 
has not previously been considered in detail. This is be- 
cause while bulges are traditionally considered to be dust 
free, like elliptical galaxies, the dust in the galaxy disc 
attenuates bulge light, particularly from the far-side. We 
can parameterize the dust attenuation curves of Fig. 1 as 
follows: 

MU gc = M b ° ulgc - h ~ Ml - c os(z)] b3 (1) 
A/ d c isc = M d ° isc - d x - d 2 [1 - cos(t)] d3 (2) 
where M° and M c represent the observed and corrected 
magnitudes, respectively, i is the disc inclination, and the 
coefficients b\ 1 62, 63 and d\, d 2 , d^ are listed in Table 1. 

Using the Millennium Galaxy Catalogue data (MG C ; 
iLiske et al.ll2"00llDriver et al.ll2005HAllen et al.ll200l we 



now derive the total B-band GLF incorporating the im- 
pact of dust attenuation. This involves seperating all 
galaxies into their bulges and disc components, correct- 
ing their fluxes according to the formulae above, and 
then rederiving their combined fluxes (note: pure ellipti- 
cal systems are not modified). Fig. 2 shows the raw ob- 
served GLF (red dotted line), the inclination corrected 
GLF (green dashed line), and the fully dust-corrected 
GLF (blue solid line, i.e., corrected for both the empiri- 
cally verified inclination-dependent attenuation plus the 
residual face-on attenuation determined by the model). 
The corresponding Schechter function fits to the GLFs, 
obtained u sing a standard step- wis e maximum likelihood 
estimator (Efstath iou et alj Il988f ). are also shown on 
Fig. 2 and tabulated in Table 2. The changes due to the 
dust correction are significant. Going from the original, 
observed -B-band GLF to the final, fully dust-corrected 
GLF, we find a 19a shift in the characteristic luminos- 
ity (M*), and a 5cr change in the faint-end slope (which 
determines the space density of dwarf galaxies) . Nearby 
galaxies therefore produce far more photons in the 400- 
450 nm range (.B-band) than previously supposed. In 
fact, after first removing the contrib ution to the GLF du e 
to ellipticals (0.14 h m L Q Mpc" 3 ; IDriver et al1l2005f) . 
only (58±5)% of -B-band photons escape from the nearby 
spiral galaxy population into the IGM (or 60 ± 5% if one 
includes the ellipticals). 

2.1. Extrapolating the impact of dust on the GLF to 
other bandpasses 

The above result provides us with a single dust cor- 
rected datum for the CSED. Unfortunately we cannot re- 
peat this measurement at other wavelengths because we 
lack the appropriate imaging data at this time. We can 
however resort to a simplification: For all possible B/T 
values from to 0.8 we derive, using the equations given 
above, the implied photon escape fraction integrated over 
cos(i). We then adopt as the effective mean B/T, that 
model galaxy whose photon escape fraction is the same 
as that determined from our GLF analysis. This yields 
an intrinsic (B/T) = 0.13±^|. Essentially, this is an 
effective average B/T value that corresponds to the vol- 
ume and luminosity weighted average of the individual 
escape fractions of all the galaxies in our sample in the 
absence of dust. 

The above has provided us with an effective average 
galaxy that represents the galaxy population at large 
in the -B-band. To now derive the equivalent effective 
average galaxy at other wavelengths we need to know 
the mean bulge and disc colours in order to modify the 
bulge-to-total ratio accordingly. For example, galaxy 
bulges are typically red (relative to the disc) so as we 
move towards longer wavelengths we expect the canon- 
ical (B/T) ratio to rise a little. To obtain the mean 
bulge and disc colours we supplement our B-band MGC 
data with multi- wavelength data pr ovided by the over- 
lappi n g Sloan Dig ital S ky Survey ^Driver et al.l 120051 . 
120071 ): (Liske et al.ll2003l I Allen et al.l l2006h . Using the 
median colours for bulge-only or disc-only systems we 
transpose our canonical bulge-to-total ratio to the SDSS 
bandpasses: (B/T) = 0.11 (u), 0.14 (g), 0.14 (r), 0.14 
(i), 0.16 (z). For JHK we adopt the 2-band (B/T) ra- 
tio and for the UV range we adopt (B/T) = (which 
implicitly assumes that star-formation has ceased in the 



Cosmic Energy Spectrum 



3 



bulge regions). 

2.2. The photon escape fraction 

We can now derive the mean photon escape fraction 
at any wavelength using the bulge and disc attenuation- 
inclination relations predicted by our calibrated dust 
model (see Section 2) coupled with the above (B /T) val- 
ues. Fig. 3 displays the corresponding photon escape 
fractions for a variety of bandpasses averaged over all 
viewing angles. To understand how critically this de- 
pends on the adopted (B/T) values we also show on 
Fig. 3 the photon escape fractions that result from as- 
suming, in each bandpass, the extreme values of (B/T) — 
(pure disc, upper dotted line) and (B/T) = 0.35 
(ear ly/mid- type disc galaxy, lower dotted line), respec- 
tively. Evidently, the photon escape fractions do not 
depend critically on the assumed "canonical" values of 
(B/T) and hence our shortcut method should be consid- 
ered robust (i.e., Fig. 3 is essentially a direct prediction 
from the MGC calibrated dust model with little depen- 
dence on (B/T)). 

2.3. The Cosmic Energy Spectrum 

The values shown in Fig. 3 can be used to derive 
the CSED corrected for dust attenuatio n. A com- 
pend iu m of recent GLF mea su rements (iDriver et ail 
20071: iBudavari et all 120051: iBlanton et all 120031: 
Kochanek et all 120011: iBell et all [20031: iBabbedee et al l 
2006L iHuang et all 120071 : iTakeuchi et all l2006f ) based 
on nearby samples (all allegedly complete, corrected to 
rcdshift zero, and converted to Hq — 70 km s _1 Mpc -1 , 
but uncorrected for dust attenuation) are shown in 
Fig. 4. Where necessary these have been converted 
from luminosity density units to energy density units 
and together they span the wavelength range from 
the far-UV to the far-IR. Also shown on Fig. 4 is our 
single fully dust-corrected B-band luminosity density 
value (from Table 2), which lies significantly above the 
previous uncorrected estimates. 

Using the photon escape fractions from Fig. 3 we now 
correct the attenuated GLF measurements to produce, 
for the first time, the unattenuated CSED (grey data 
points on Fig. 4). This constitutes the actual spectral 
energy output in the Universe today due to the total inte- 
grated starlight before dust attenuation. To calculate the 
total energy of starlight versus that which escapes into 
the IGM we need to integrate over these two datasets. In 
order to interpolate across the full far-UV, optical, and 
IR wavelength range we adopt a recent stellar s ynthe - 
sis model (PEGACE, [F ioc fcR occa-Volmeragenl (|1997f ) 
see iBaldrv fc Glazebroo k ( 2003) for details of the mod- 
elling) which provides a reasonable fit to both the attenu- 
ated (orange line) and, when corrected, the unattenuated 
(black line) CSED. Integrating these two curves yields a 
total stellar energy output of (1.6 ± 0.2) x 10 35 W Mpc -3 
of which (0.9 ± 0.1) x 10 35 W Mpc -3 escapes into the 
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IGM. Note that these numbers take into account a 10% 
compone nt of the CSED lq ngwards of 400 nm due to el- 
lipticals (jDriver et al.ll2005l ). For the CSED to be in equi- 
librium (and energy conserving) the difference between 
these two energy values, (0.7 ± 0.2) x 10 35 W Mpc -3 , 
should not exceed the total emission from dust in the 
far-IR, the traditional sticking point for optically thin 
models. Fig. 4 also shows a model of the dust emission 
(red line; iDale fc Heloull2002T ) that reproduces fairly well 
the observed data Integrating the far-IR curve we obtain 
a total radiant dust energy of (0.6± 0.1) x 10 35 W Mpc -3 , 
in excellent agreement with our prediction. This provides 
independent support that the high value for the central 
face-on opacity is correct, and implies that significant 
corrections to the observed flux of individual galaxies 
are necessary. This result also leaves little or no room 
for other sources of dust heating in the nearby Universe, 
such as active galactic nuclei, and provides the first fully 
reconciled estimate of the CSED. 

3. DISCUSSIONS AND RAMIFICATIONS 

This work has reconciled three apparently inconsistent 
observations, namely; the severe at tenuation- inclinatio n 
relation seen in the MGC data (jDriver et all l200l : 
the conclusion tha t inter-arm regions are optically thin 
(| White et all 120001: iHolwerda et al.| [2007): and the rela- 
tively high far-IR dust emission. The TP-model achieves 
this by incorporating distinct dust components that al- 
low for optically thin inter-arm and outer regions coupled 
with an optically thick central region. This relatively 
simple development, has a number of far-reaching rami- 
fications. Firstly, all basic measurements of galaxy fluxes 
that do not correct for dust attenuation will require sig- 
nificant revision (i.e., 0.2-2.5 mag in B) depending on an 
individual galaxy's inclination, bulge-to-total ratio and 
wavelength of observation. Second, many galaxies will 
contain heavily embedded bulges due to the centrally 
concentrated dust in their discs. This can easily lead 
to significant errors in optical estimates of their fluxes 
and stellar masses; this is because although stellar mass 
estimates do correct for optically thin dust attenuation 
they cannot, of course, correct for mass hidden behind an 
entirely optically thick screen. Finally, dust attenuation 
could conceivably play a part in the morphology-density 
relation and the proposed transformation of disc galaxies 
from late- to early-type as they enter the cluster environ- 
ment. For example, using the corrections provided, it is 
easy to show that an Sab galaxy at the median inclina- 
tion (60°) will see its observed B/T change from 0.3 to 
0.6, its color get significantly redder, and its luminosity 
get brighter, if all its dust was to be removed. 
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TABLE 1 

Coefficients for use in Equations 1 and 2 to provide bulge and disc attenuation corrections (with < 5% error in the 

predicted values due to the adopted fitting function). 



Bandpass 


bi 


f>2 


&3 


di 


d 2 




u 


1.10 


0.95 


2.18 


0.45 


2.31 


3.42 


B 


0.89 


1.27 


1.73 


0.24 


1.20 


2.73 


9 


0.83 


1.29 


1.71 


0.22 


1.18 


2.74 


r 


0.63 


1.33 


1.73 


0.16 


1.10 


2.80 


i 


0.48 


1.35 


1.84 


0.11 


1.03 


2.89 


z 


0.38 


1.35 


1.84 


0.09 


0.96 


2.98 


J 


0.25 


1.22 


2.26 


0.06 


0.80 


3.21 


H 


0.18 


1.02 


2.43 


0.05 


0.64 


3.51 


K 


0.11 


0.79 


2.77 


0.04 


0.46 


4.23 



TABLE 2 

Derived Schechter luminosity function parameters for the MGC with varying degrees of dust attenuation corrections. 



Sample 


M* — 5 log h 7(] 


a 


4>* 


is 




[mag] 




[10~ 3 Mpc" 3 (0.5 mag)" 1 ] 


[10 8 h m L Q Mpc- 3 ] 


No dust corr. 


-20.57 ±0.04 


-1.14 ±0.03 


6.7 ±0.3 


1.9 ±0.2 


Inclination corr. 


-20.78 ±0.04 


-1.16 ±0.03 


7.0 ±0.3 


2.4 ±0.2 


Incl.+face-on corr. 


-21.32 ±0.05 


-1.32 ±0.02 


4.8 ±0.3 


3.1 ±0.5 
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Fig. 1. — The dust attenuation— inclination relations for galaxy discs (right) and bulge s (left) shown for a variety of bands (as indicated 
on the right hand sides) as predicted by the TP-model calibrated on B-band disc data (Driver ct al. 2007]) • The non-zero attenuation at 
1 — cos(i) = is the residual face-on attenuation predicted by the model. 
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Fig. 2. — Main panel: the B-band galaxy luminosity function; ignoring all consideration of dust (dotted line), after consideration of the 
empirically derived attenuation— inclination relation only (dashed line), and after a full treatment of dust attenuation (including the face-on 
model correction; solid line). Side panel: the projected 3<r error contours for two of the three fitting parameters, a (the faint-end slope) 
and M* (the characteristic luminosity). 
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Fig. 3. — The global photon escape fraction averaged over all inclinations versus wavelength. The dotted lines show estimated error 
boundaries (see text). 
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Fig. 4. — The cosmic energy output covering the region dominated by starlight (left peak) and by dust emission (right peak). The orange 
line shows the observed (uncorrected) cosmic energy output from the total nearby galaxy population, while the black line shows the same 
after correction for the fraction of photons attenuated by dust. The discrepancy in the integrals over these two curves yields the total 
energy of starlight lost to heating of the dust grains. If starlight is the only source of dust heating then this energy loss must equal the 
total radiant energy of the dust emission (i.e., the two shaded regions must and do contain equal energy). 



